CLAIMS 

1.-52. (CANCELLED) 

53. (CURRENTLY AMENDED) A method of estimating arterial delay and arterial 
dispersion (t, a, a) values for outputting blood perfusion indices for a region of 
interest (ROI) by #em operating a computer program on intensity data [[in]] input 
to a computer comprising: 

a. using a computer to apply applying a first gamma-yariate function 
(GVF) to an arterial input function (AIFa) using a comput e r to provide an 
estimated first model of a vascular transport function ha(t), wherein for t < 

ti, ha(t) = 0 and for t > ti, ha(t) = — - t^Y' ^'O'hV^i ^ ^y^^^^^ is the 

CTi 

transit time of a contrast agent from a measured initial said AIFa to a 
region of interest (ROI) and 0\ is an estimatin g estimate of an initial d e lay 
dispersion value of said contrast agent, wherein said Oi=(ti)(Pi)/(l-Pi), 
wherein said Pi is a known relative dispersion value having a range from 0 
to 1; 

b. using a comput e r to convolv e convolving AIFa(t) with said ha(t) with 
a 1=0 using a comput e r for obtaining an arterial input function AIFt(t) = 
AIFa(t) ® ha(t) with ai=0 at said ROI; 

c. using a computer to e stimat e estimating a blood flow rate Ft and a tissue 
impulse residue function Re(t) using a computer by deconvolving a 
concentration curve C(t) = (Ft/kH)AIFt(t) ® Rc(t), wherein kn is a 
hormatooH t hematocrit correction constant having a known value , wherein 
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a tissue transport function hp(t) is determined by tuft) = - dIL(tVdt: 

4: using a comput e r to optimize said m e an transit tim e and disp e rsion (t^ ? 

e^-O a) valu e s using a l e ast squar es m e thod from an e stimated transport 
function b^(t); and 

e;^ using a computer to output inputting said estimated and optimiz e d 

calculat e d tissue m e an transit tim e and dispersion fti, ei b zjefa ) valu e s from 
an e stimated transport function he(t) for input to a simulated transport 
function h^(t), wh e r e in a simulat e d tissue impuls e r e sidue function R^ (^4s 
d e t e rmined, wher e in a simulat e d conc e ntration curve Cs(t) is fitted to said 
measur e d C(t) and quantitativ e said blood perfusion indic e s are calculated, 
wh e r e in each said stop is p e rform e d by a suitably programm e d computer. 

d. determining a simulated transport function hj[t)= — (r-^)"'^'^^"^'^^''^ 



when t^t? and h.(t')=0 when t<t2, wherein A. = a2***T(l + an) , wherein said 

h^(t) = — when said t2=0. or said hjt) = — g-<^-^2)/^2 ^^en said 
At 02 

g? = 0. wherein said o^^is a dispersion of said hjt) and said t? is a mean 
transit time of said h^(t), wherein when said a? = 0 a peak height 
fPH>l/a2_and a mean transit time fMTTl=t?+a? are used to determine 
said g^: 

e. determining a simulated tissue impulse residual function (IRF) 



R^(t) = 1 - J h/r)dr . wherein a simulated contrast agent concentration 



0 



am is determined bv CJt) = (F^ lk^)AIFXt) ® RJt) : and 
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f. fitting said simulated Cffl to said C(t) by iteratively minimizing S 
using a least squares method defined by 5 = ^(Cit) - C/t))^ . wherein 

said iteratively mimimizing S comprises reducing the number of 
adjustable parameters, wherein said adjustable parameters are reduced to 
5 fiye by fixing said a i ^=0 and said h = 0, or by fixing said aj^=Q and said 

CC2=0. and wherein said adjustable parameters are further reduced to four 
by fixing said relative dispersion Pi= 0\/(0}H}) of sai d hJt) resulting in 
said oi dependent on said ti; 
wherein each said step is performed by a suitably programmed computer. 

10 

54. (CURRENTLY AMENDED) The method of claim 53, wherein said intensity 
data is generated by administering a contrast agent to a body lumen of a body 
during a dynamic imaging scan, wherein said body lumen comprises an artery 
or a_yein, wherein an image response fi-om sai^ contrast agent is recorded to 
1 5 computer data storage in a computer. 

55. (PREVIOUSLY PRESENTED) The method of claim 53, wherein said C(t) is a 
temporal concentration of said contrast agent obtained from said intensity data, 
wherein said intensity data comprises contrast images sequentially acquired 

20 from a region in a body, whereby said contrast agent concentration is plotted 

versus time. 

56. (CURRENTLY AMENDED) The method of claim 53, wherein said AIFa is 
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based on a measured e arly arrival contrast ag e nt p e ak intensity profile from a 
feeding blood vessel to said ROL 

57. (CURRENTLY AMENDED) The method of claim 53, wherein said AIFa is 
5 scaled upward according to a venous input function (VIF), wherein said VIF is 

based on a measured lat e arrival contrast agent peak intensity profile from a 
large vein draining firom said ROL 

58. (PREVIOUSLY PRESENTED) The method of claim 53, wherein said 
10 estimated transit time ti is the transit time of said contrast agent fi^om a 

measured mitial said AIFa of said contrast agent C(t) in a body lumen to said 
ROI, wherein said ti is estimated from plots of said AIFa versus time and said 
C(t) versus time. 

15 59. (PREVIOUSLY PRESENTED) The method of claim 53, wherein said ha(t) is 

calculated using said estimated transit time ti and said estimated dispersion 
value Oi, wherein ha(t) with ai=0 is plotted versus time. 

60. -65. (CANCELED) 

20 

66. (CUflRENTLY AMENDED) The method of claim 53, wherein said AIFt(t) is 
measureable in a s mall lumen showing a delay relative to said AIFa(t), wherein 
optimized values for said 0\ and said ti are determined by fitting said ,simulatod 
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convolved A IFt(t) to said measured AIFt(t), wherein said relative dispersion Pi 
is determined and applied to all other said intensity data of said ROI using said 
Pi, wherein a more robust fitting process is provided by a reduced number of 
parameters for optimization. 

67. (CURRENTLY AMENDED) The method of claim 66, wherein when said 
relative dispersion Pi is determined, said vascular transport function ha(t) is 
described by a single variable said ti with a constant said Pi, wherein a 
two-step method is used to determine said delay and said dispersion values 
comprising: 

a. deriving an initial tissue impulse residue function Ro(t) by 
deconvolving C(t) = (Fo/kH)AIFa(t) ® Ro(t) using a model-free singular 
valu e d e composition (SVD) deconvolution m ethod, wherein said time 
delay ti is determined by a maximum position of said Ro(t) at Ro 
max(t=ti); and 

b. determine said AlFt(t) at an input of said ROI using said ha(t) with said 
ti and said Pi held constant, wherein said Ci is determined. 

68. (CURRENTLY AMENDED) The method of claim 67, wherein a value 
of tissue blood flow Ft and a corrected impulse residue function Re(t) 
are obtained by deconvolving C(t) = (Ft/kH)AlFt(t) ® Rc(t) using said 
SV D model-free deconvolution method, wherein said perfusion indices 

poo 

are determined from a curve of said Rc(t), wherein MTT= J Re(T)dT, 
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BF=Ft, andBV=BF*MTT. 



69. (CURRENTLY AMENDED) The method of claim 53, wherein said contrast 
agent is in a tissue ROI having a tissue mean transit time t, wherein a tissue 
impulse residue function is approximated by the relation R(t >t) = Ee'^^^'"^ and 
R(t^T) = 1, wherein E is an extraction fraction of said contrast agent in said 
tissue, wherein k is a constant clearance rate of said contrast agent diffusing 
from said tissue having a relation k = E*FtA^e, wherein Ve_is the volume 
fraction of extravascular and extracellular space (EES) in said tissue, 

70. (CURRENTLY AMENDED) The method of claim 69, wherein said tissue 
impulse residue function Rs(t) of said simulated concentration curve Cs(t) is 
replaced by an average impulse residue function that incorporates said 
contrast agent leaked out of a blood vessel into said tissue and gradually 
clearing from said tissue, wherein said simulated concentration curve Cs(t) 
is fitted to said measured C(t) and quantitative said blood perfiision indices 
are calculated, wherein said E and said Ve are additional parameters 
optimized with other adjustable parameters using a least squares method. 
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